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摘要 ”本 文通 过 对 历史 太阳 周 的 峰值 、 上 升 期 、 下 降 期 等 特征 参数 的 统计 分 析 ,， 给 出 了 峰值 变化 与 上 
升 期 之 间 的 三 类 线性 拟 合 关 系 ， 两 个 太阳 周 峰 值 之 间 长 度 与 下 降 期 的 两 类 拟 合 关系 ,并 通过 分 析 太 阳 周 
的 长 期 变化 周期 和 第 10 周 之 后 的 偶 、 奇 数 峰值 的 变化 规律 ,给 出 第 25 太阳 周 的 预测 结果 : 预计 第 25 太 
阳 周 在 2020 年 6 月 前 后 开始 , 该 周 会 比 第 24 太阳 周强 , 太阳 黑子 数 峰值 为 119.2+5.5, 峰值 强度 类 似 于 
第 17 太阳 周 和 23 太阳 周 . 


关键 词 。 太阳 周 , 峰值 ， 上 升 期 ， 下 降 期 
PACS: 96.50.Wx, 96.60.-j, 96.60.Q-, 96.60.qd 


j 
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1 引言 经 常用 于 表示 太阳 活动 水 平 的 指数 之 一 , 它 是 由 日 
掉 上 的 黑子 群 数 与 黑子 个 数组 合 得 到 凹 ,尽管 它 不 如 


空间 环境 对 航天 活动 的 影响 越 来 越 受到 技术 部 。 太阳 在 2800 MHZ 的 辐射 流量 密度 或 时 子 群 面积 总 
门 的 关注 , 无 论 是 航天 器 设计 、 发 射 和 回收 , 还 是 飞 数 思 等 指数 那样 具有 明确 的 直观 性 和 物理 意义 但 它 


行 计 划 编 制 、 轨 道 飞 行 姿 态 控制 操作 等 ， 都 需要 考虑 有 长 达 二 百 多 年 的 观测 历史 ,作为 一 个 统计 量 能 清 
空间 环境 的 影响 . 太阳 活动 是 空间 环境 的 扰动 源头 ， 楚 地 反映 出 太阳 活动 的 周期 性 外 ,作为 太阳 活动 区 的 
是 引起 空间 环境 灾害 性 事件 的 主要 原因 . 太阳 活动 ” 核心 成 员 它 与 太阳 次 斑 申 日 黑 物 质 抛射 9 等 活动 现 
水 平 的 高 低 直 接 影响 到 近 地 空间 环境 状态 (大 气 环 象 的 产生 有 较 好 的 统计 相关 ， 因而 它 比 其 他 太阳 指 
境 、 电 离 层 环境 等 ). 在 太阳 活动 剧烈 变化 的 情况 下 ， ” 数 更 多 地 应 用 于 太阳 活动 物理 、 空间 天 气 和 地 球 物理 
大 气 密度 可 达到 数 倍 的 变化 幅度 ,严重 影响 正常 的 ”等 领域 . 在 分 析 太 阳 黑 子 数 的 长 期 变化 、 评 价 长 期 太 
轨道 维持 . 因此 对 太阳 活动 的 准确 预报 是 空间 环境 阳 活 动 水 平时 , 广泛 使 用 的 是 太阳 黑子 数 月 均值 的 
预报 的 重要 课题 之 一 . 平滑 值 ( 即 平滑 月 均值 ). 该 值 是 国际 上 界定 太阳 周 包 

在 太阳 活动 预报 研究 中 , 沃 尔 夫 太 阳 黑 子 数 是 ”位 相 和 黑子 数 极 值 的 参量 , 因此 也 是 太阳 周 预报 的 


bhol 
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由 于 太阳 黑子 数 可 能 存在 约 有 25% 的 不 确定 
度 中 ,对 太阳 黑子 数 的 预报 仍然 是 一 个 十 分 艰巨 的 任 


务 . 近 百 年 来 研究 者 采用 


了 各 种 各 样 的 方法 对 太阳 
活动 (特别 是 峰值 ) 进 行 预报 : 黑子 数 长 短 周期 法 、 黑 
子 数 峰值 变化 趋势 法 、 奇 / 偶 周 法 、 时 间 序 列 法 、 地 
磁 先 兆 法 、 相 似 周 法 及 发 电机 模型 等 . 在 第 23 太阳 


周 以 前 ,由 于 先兆 方法 在 连续 几 周 (20-22) 进 行 了 较 
为 准确 的 预报 ， 因 而 这 种 方法 被 普遍 接受 ,由 此 也 促 


进 了 太阳 发 电机 模型 的 发 


展 "但 这 种 方法 在 第 23 


周 却 失效 了 , 得 到 的 预报 结果 普遍 偏 高 . 科学 工作 者 


为 了 提高 预报 效果 提出 了 


些 新 的 参量 ， 如 拐点 处 


的 斜率 中、 黑子 相 对 数 与 地 磁 指 数 的 相对 变化 率 
a Li) Kane02 和 Pishkalo 等 人 03 曾 对 第 22 周 、23 


周 及 24 周 的 预报 结果 进行 了 比较 , 可 以 看 到 ， 由 于 


式 中 Mi 为 由 观测 得 到 的 第 i 月 的 太阳 黑子 数 月 均值 . 
一 个 太阳 周 中 ，M, 的 极 小 值 及 极 大 值 常 分 别 记 为 


Main 和 M nax 


， 上 升 期 的 时 间 长 度 记 为 R， 下 降 期 的 


时 间 长 度 记 为 D. 这 些 量 和 M, 的 时 变 曲线 的 形态 特 


点 ， 共 同 描述 了 一 个 太阳 周 的 太阳 活动 水 平 的 演化 
概貌 . 为 方便 统计 ， 本 文中 定义 了 相 邻 两 个 太阳 周 分 
别 为 第 nn AME ntl 周 的 特征 参数 ， 如 图 1 所 示 . 第 


n 周 的 峰值 为 MM;,(n)， 上升 期 为 RCn)， FEN Din), 


第 n+l 周 的 峰值 为 Mi (ntl), LAA Rent), 两 


个 太阳 周 峰 值 M an), ML, (+1) TE EAD IY Ti) 4G FE 
示 为 K(n, n+1). 
本 文 统计 过 程 中 采用 的 数据 是 国际 公开 发 布 的 


太阳 周 黑子 数 相关 特征 参数 . 国际 上 沃 尔 夫 太阳 黑 


研究 者 采用 方法 的 不 同 ,对 同一 参数 的 预报 结果 差 


异 也 比较 大 ， 比 如 第 24 周 的 太阳 黑子 数 峰 值 预报 结 
果 最 小 为 50, 最 大 达到 190. 通过 大 量 的 预报 方法 和 


预测 效果 比 对 来 看 ， 目 前 还 没有 一 种 方法 对 最 近 的 
儿 个 太阳 周 都 做 出 了 较 准 确 的 预报 ， 太阳 周 峰值 预 


报 仍然 是 一 个 有 待 解 决 的 问题 . 


通过 统计 分 析 以 往 各 太阳 周 参数 间 的 变化 特征 ， 
从 中 寻找 出 与 未 来 太阳 周 相 关 的 先兆 因子 并 进行 预 


测 ， 这 在 太阳 周 长 期 预报 方法 中 最 为 常见 . Wang04， 
Lil), Du 等 人 09 通 过 统计 分 析 上 升 期 与 峰值 的 关系 


预测 第 24 周 的 峰值 及 时 间 ; 


Jain''”, Hathaway 等 人 103 


分 别 开 展 了 地 磁 与 太阳 黑子 的 变化 统计 分 析 ; Du 等 
人 通过 统计 分 析 提 出 了 一 种 加 权 方 法 预测 峰值 ， 
Badalyan ”开展 冕 绿 线 与 太阳 黑子 的 相关 性 来 进行 


黑子 数 预 测 . 本 文 对 典型 太阳 周 特征 参数 间 的 变化 


关系 进行 统计 分 析 , 并 按 其 变化 特征 进行 不 同 分 类 
和 公式 拟 合 ， 最 后 给 出 第 25 周 预报 结果 ,并 与 已 发 


表 的 预报 结果 进行 比较 . 


2 数据 和 统计 方法 


对 一 个 太阳 周 中 每 个 月 的 太阳 活动 水 平 的 预报 ， 


是 以 太阳 黑子 数 平滑 月 均值 作为 预报 量 . 式 (1) 给 出 


了 第 1 月 黑子 数 平滑 月 均值 于 ,的 定义 
1 1 i+5 
AET pd at MM, 2 () 


j=i-5 


max 


子 数 月 均值 从 1749 年 1 月 开始 记录 , 根据 公式 (1)， 


太阳 黑子 数 第 一 个 平滑 月 均值 从 1749 年 7 月 开始 记录 ， 


而 国际 上 第 一 个 太阳 周 是 从 1755 年 3 月 开始 算 起 ， 目 


前 处 于 第 24 太阳 周 . 本 文 统计 数据 从 1749 年 7 月 起 ， 


为 便于 说 明 ， 


Wang 等 人 ! 


定义 该 周 为 第 0 太阳 周 . 在 2002 年 ， 
4 曾经 统计 过 峰值 与 上 升 期 之 间 的 关系 ， 


并 按照 上 升 速度 划分 为 高 上 升 速度 和 低 上 升 速度 两 类 ， 


统计 中 不 包含 第 0 太阳 周 . 本 文通 过 对 0-23 周 上 升 
期 Rn) 与 峰值 M,,(n) 的 统计 关系 发 现 ，24 个 周 的 峰 
值 与 上 升 时 间 的 分 布 特征 可 以 通过 三 类 拟 合 关系 来 反 
Bk: 第 一 类 包括 1, 5, 7, 19 Al, 第 二 类 包括 0, 9, 21 周 ， 
其 余 太 阳 周 属于 第 三 类 . 图 2 给 出 了 三 类 线性 


太阳 黑子 平滑 月 均值 


RN “DM Rn 


M 


max 


(n+1) 


K(n, n+1) ——+: 
时 间 (a) 


1 太阳 周 特征 参数 定义 


Figure 1 Definitions of parameters for solar cycle. 
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2.0 2.4 2.8 3.2 3.6 4.0 4.4 4.8 5.2 5.6 6.0 6.4 6.8 7.2 
R(n) (a) 


2 (网 络 版 彩 图 )0-23 太阳 周 上 升 期 R(n) 与 峰值 Mnn) 
的 变化 关系 


Figure 2 (Color online) Relationship between R(n) and Mpa (n) 
for n=0, 1, 2, =, 23. 


拟 合 关系 及 相关 系数 ， 可 以 看 到 , 三 种 线性 拟 合 都 具 
有 很 好 的 相关 性 ,相关 系数 分 别 为 -0.998，-0.984 和 
-0.963. 按照 三 类 线性 拟 合 关系 ,本 文 给 出 了 0-23 周 
峰值 Ma (M 随 上 升 期 R(n) 的 线性 拟 合 公式 . 
第 1, 5,7 和 19 周 峰值 与 上 升 期 变化 符合 拟 合 公 
式 (2), 拟 合 结果 与 实测 值 的 均 方 根 差 为 AM=1.46. 

M „a (n) = 362.65 — 44.53R(n). (2) 
第 0, 9 和 21 周 峰 值 与 上 升 期 符合 拟 合 公式 (3)， 
拟 合 结果 与 实测 值 的 均 方 根 差 为 AM=6.11. 
M „ax (n) = 317.24 — 42.47R(n). (3) 
其 余 太 阳 周 (包括 第 2-4, 6, 8, 10-18, 20, 22 Fil 23) 
变化 符合 拟 合 公 式 (4), 拟 合 结果 与 实测 值 的 均 方 根 
差 为 AM=5.46. 


Ms (1) = 263.33—37.92R(n). (4) 


通过 上 面 统计 分 析 ， 本 文 将 目前 已 有 的 0-23 个 
太阳 周 按照 峰值 与 上 升 期 的 关系 划分 为 三 种 类 型 的 
变化 并 给 出 相应 拟 合 公式 . 如 知道 一 个 太阳 周 上 升 
期 ， 便 可 以 通过 相关 拟 合 公式 预测 其 峰值 的 变化 范 
围 . 根据 公式 (2)-(4)， 要 预测 第 25 周 峰值 ， 就 需要 知 
道 第 25 周 上 升 期 长 度 . 为 确定 第 25 周 的 上 升 期 , 本 
文 对 R(n), D(n) 及 K(n, n+1) 各 个 参数 进行 了 统计 ,图 
3 给 出 了 0-23 周 中 K(n, n+1) 随 D(n) 分 布 关系 . 可 以 
看 到 下 降 期 Dn) 和 Kin, n+1) 具 有 很 好 的 相关 性 ， 并 
通过 线性 拟 合 关系 发 现 ， 用 两 类 拟 合 关 系 来 表示 其 


K(4,5) 


17L * CC=0.896 
[ CC=0.995 


K(n,n+1) (a) 


3.5 40 45 50 55 60 65 70 75 30 35 90 9.5 10.010.51 1 0 
D(n) (a) 
3 (网 络 版 彩 图 )0-23 太阳 周 下 降 期 D(n) 与 K(n, n+1) 的 


变化 关系 


Figure 3 (Color online) Relationship between D(n) and K(n, n+1) 
for n=0, 1, 2, …, 23. 


分 布 更 为 合适 . 第 一 类 为 K(0,1), K(4,5), K(5,6) F 
K(6,7)， 相 关系 数 为 0.995， 其余 属于 第 二 类 ,相关 系 
BA 0.896. 按照 两 类 线性 拟 合 关系 ， 本 文 给 出 了 
0-23 周 下 降 期 D(n) 的 与 Kn, n+1) 线 性 拟 合 公 式 . 

当 n 为 0,4,5 和 6 时 , D(n) 与 K(n,n+1) 的 变化 符 
合 拟 合 公式 (3)， 拟 合 结果 与 实测 值 的 均 方 根 差 为 
AT=1.11 (a). 


K(n,n+1) = 4.09 +1.23D(n). (5) 

“4 nH 1-3, 7-23 时 , D(n) 与 Kln,n+1) 的 变化 符合 

拟 合 公 式 (60),， 拟 合 结果 与 实测 值 的 均 方 根 差 为 
AT=1.03 (a). 


K(n,n+1) =2.3341.24D(n). (6) 


3 第 25 太阳 周 预测 结果 和 讨论 


通过 统计 分 析 , 第 二 节 中 给 出 了 第 n 活 动 周 峰 值 
M,,(n) 随 上 升 期 R(n) 的 拟 合 关 系 , 第 n 周 下 降 期 
D(n) 与 第 n, n+1 周 峰 值 间 时 间 长 度 K(n, n+1) 的 拟 合 
KA. 即 知 道 了 第 n 周 的 下 降 期 , 便 可 以 预测 第 n+l 
周 的 上 升 期 (Rn+1)=K(n, n+1)-D(n))， 从 而 可 以 预测 
n+1 周 的 峰值 大 小 . 

第 24 太阳 周 开 始 于 2008 年 的 12 月 ,截至 2014 
年 12 月 , 已 经 过 了 6 年 ， 从 已 发 展 的 态势 来 看 ， 第 24 
周 的 确 是 一 个 比较 弱 的 太阳 周 , 太阳 黑子 平滑 月 均 
值 在 2012 年 2 月 达到 一 个 小 峰值 66.9 之 后 就 开始 现 
下 降 , 之 后 近 1 年 都 在 58 附近 振荡 ， 从 2013 年 5 月 
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黑子 数 平滑 月 均值 预报 结果 , 81. 


开始 缓慢 上 升 ， 2014 年 4 月 达到 第 24 太阳 周 以 来 的 
最 高 值 81.9, 之 后 的 5-8 月 逐渐 下 降 ， 
79.7, 78.6 和 75.6. 那么 第 24 周 的 峰值 是 81.99 根据 
前 国际 、 国 内 预报 机 构 (NASA, IPS，SEPC) 的 太阳 
各 是 第 24 周 的 峰值 ， 


9} 


分 别 为 80.5, 


之 后 将 进入 第 24 周 下 降 期 . 图 4 是 我 们 利用 相似 周 


TRE? HM IES 24 周 的 平滑 月 
测 结果 比较 ， 从 图 中 可 以 看 到 ， 


开始 振荡 下 降 . 


均值 的 变化 剖 线 与 观 
我 们 预测 的 峰值 小 于 
实测 结果 ,峰值 出 现时 间 为 2014 年 5 月 , 之 后 黑子 数 


按照 第 24 周 峰值 为 81.9 的 情况 ,第 24 周 的 上 升 
时 间 为 5.4 年 , 那么 按照 峰值 随 上 升 期 的 分 布 特征 ， 


第 24 周 更 加 符合 图 2 中 的 第 二 类 线 怕 


E 拟 合 关 系 ， 按 


照 对 应 拟 合 公 式 (3), 计算 峰值 为 87.90+6.11, 与 已 经 


出 现 的 81.9 相近 . 


对 于 第 24 周 的 下 降 期 预测 , NASA F IPS 的 预测 


结果 相近 ,NASA 预测 在 2019 4 


在 2019 年 12 月 结束 . 本 文 按照 相似 周 特征 参数 相似 
原则 R21 取 第 24 周 的 相似 周 (第 12, 14 周 ) 的 平均 下 降 
期 为 第 24 周 的 下 降 长 度 ， 即 为 (6.2+0.4) 年 ,由 第 二 节 
中 D(n) 与 K(n, n+1) 的 两 类 线性 拟 合 关 系 可 以 看 到 ， 
只 有 4 个 周 属于 第 一 类 ， 且 都 处 在 第 7 周 之 前 , 其 后 
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结束 , IPS 预测 


第 7-23 周 均 符 合 第 二 类 线性 拟 合 特 征 . 因此 ， 本 文 
WA n N 24 时 K(24, 25) 的 变化 应 符合 第 二 类 拟 合 关 
A. 因此 应 用 公式 (5) 可 得 到 K(24, 25) 为 10.02+0.5. 
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4 (网络 版 彩 图 ) 相 似 周 方法 预测 的 第 24 周 黑子 数 与 观 


测 值 比较 


Figure 4 (Color online) Comparison of observation and prediction 


for cycle 24. 


再 由 K(24, 25)-D(24) 可 预测 第 25 周 上 升 期 为 3.8+0.1， 
即 第 25 周 预计 起 始 于 2020 年 6 月 前 后 ,峰值 将 在 
2024 年 4 月 前 后 出 现 . 
得 到 上 升 时 间 ， 即 可 通过 上 升 期 和 峰值 的 拟 合 
公式 (2)-(4) 对 应 计算 得 到 第 25 周 的 峰值 分 别 为 
193.4+1.5, 155.8+6.1 和 119.2+5.5. 那么 , 第 25 AY 
峰值 属于 哪 类 预报 结果 呢 ? 太 阳 活 动 除 了 准 11 年 的 
周期 变化 外 , 还 具有 长 周期 变化 特征 1， 这 已 获 等 
众多 研究 者 的 公认 , 且 一 些 人 认为 正 是 因为 太阳 黑 
子 受 到 了 长 周期 的 调制 , 才 导 致 第 23 周 众 多 预报 结 
果 不 准 . 2010 年 HathawayP8l 拟 合 给 出 了 峰值 随 活 动 
周 数 变化 的 正弦 函数 曲线 ,并 认为 峰值 具有 9.1 个 太 
阳 周 变化 长 周期 . 而 最 近 一 个 峰值 出 现在 第 19 周 ， 
第 25 周 应 处 于 9.1 个 太阳 周 变化 长 周期 的 第 6 周 , 将 
不 会 达到 强 活动 周 水 平 , 但 会 强 于 第 24 周 . 根据 偶 
数 周 与 其 后 奇数 周 峰 值 变化 情况 (图 9)， 可 以 发 现 ， 
从 第 10 周 以 来 除了 第 23 周 呈 现 异常 之 外 ， 紧 接 偶数 
周 之 后 的 奇数 周 的 峰值 总 高 于 其 前 面 的 偶数 周 黑子 数 ， 
第 11, 13, 15, 17, 19, 21 周 缘 如 此 ,这 也 从 男 一 个 方面 
反映 了 按 规律 变化 第 25 周 峰 值 会 高 于 第 24 周 . 

综 上 所 述 ,通过 太阳 周 的 长 周期 变化 特征 、 偶 数 
周 和 奇 数 周 峰 值 的 变化 规律 及 拟 合计 算 结果 ， 本 文 
认为 , 第 25 周 不 会 太 强 , 但 峰值 要 高 于 第 24 周 , 预 
计 太 阳 黑 子 数 峰 值 为 119.4+5.5. 表 1 给 出 了 其 他 作 
g 113:19.20.29- S51 25 周 的 预测 结果 ,可 以 看 到 ,预测 结 
果 并 不 一 致 , Javaraiah22 和 Abdusamatov 中 认为 第 25 周 
峰值 在 50 左右 ，Hathaway60 预 测 在 70 £A, Dul’? 
给 出 了 102, 111 和 144 三 个 结果 . 本 文 的 预测 结果 与 
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Figure 5 Variation of maximum from Even cycle to odd cycle. 
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表 1 第 25 太阳 周 预报 结果 的 预报 方法 和 相关 因子 . 由 于 目前 的 预报 方法 主要 
Table} Predictions of solar oy 25 建立 于 对 历史 数据 的 统计 分 析 基 础 上 ,预报 方法 的 
黑子 3 | LH SUL ERY [E Z = Bie Vp Sd lip = 
— - = = 不 同 , 结果 会 有 很 大 差异 ,目前 还 没有 一 种 方法 对 最 
avaraiah, 2014 50+10 一 
Abdusamatov, 2007”! 50415 = 近 的 几 个 太 BH Jal 都 做 出 了 较 准 确 的 预报 > X 阳 周 峰 
Hathaway, Wilson, 2004°" 70+30 2023 和 直 预 报 仍然 是 一 个 有 待 解决 的 问题 . 

Bd or - 在 工程 应 用 中 ,为 了 满足 任务 需求 , 往往 需要 提 
eee 和 ADT 供 未 来 几 十 年 的 太阳 活动 预报 , 即 在 第 24 周 还 在 进 
ra 117 2024 行 中 时 ， 就 需要 知道 第 25 周 的 太阳 活动 情况 . 基于 
121 20293 此 需求 ,， 本文 对 0-23 太阳 周 的 峰值 、 上 升 期 、 下 降 

Du, 2006" 111.6417.4 - i. S a apa S 
Hamid and Galal, 2013" 118 2020-2030 分 析 . 按照 分 布 特征 ,把 峰值 随 上 升 时 间 的 关系 划分 
Rigozo, 2011 132:1 E 为 三 类 ,两 周 峰 值 之 间 的 时 间 随 下 降 期 的 关系 划分 

Du, 2006"°! 144.3427.6 一 ‘ Me eb LIA JN ye ie 
为 两 类 ， 并 给 出 相应 的 线性 拟 合 公式 . 根据 目前 第 


Quassim*”!, Hamid 等 人 中 1 的 结果 一 致 . 


24 周 的 发 展 趋势 和 预报 结果 ,以 及 太阳 周 的 长 周 基 
变化 特征 和 第 10 周 之 后 的 偶 、 奇 数 峰值 的 变化 规 得 
pa 给 出 第 25 周 的 预报 : 预计 第 25 周 在 2020 年 6 月 前 
结论 后 开始 ， 该 周 会 比 第 24 周强 ,峰值 预计 会 在 2024 年 
几 十 年 来 ,为 准确 预测 太阳 周 的 特征 参量 , 特别 4 月 前 后 出 现 , 黑子 数 为 119.2+5.5， 峰 值 强度 类 似 于 
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是 太阳 黑子 数 峰值 ,太阳 活动 研究 者 不 断 在 寻找 新 第 17 周 和 23 周 . 
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The prediction of maximum amplitude of solar cycle 25 
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In this paper, the correlation between various parameters of solar cycles 0-23 is investigated. The parameters include 
maximum, rising time and descending time. Three derived regression equations are given between maximum and 
rising time. Two derived regression equations are given between descending time and total of descending time and 
next rising time. Additionally, according to the cycle periodicity of secular trend and the regularity of maximum 
between even cycles and following odd cycles, the prediction of cycle 25 is obtained. It is expected that solar cycle 25 
will begin in June 2020 and will be relatively stronger than cycle 24. The sunspot number maximum of cycle 25 is 
predicted to occur in April 2024 with a sunspot number of 119.2+5.5, which will be nearly the same strength of cycle 
17 and 23. 


solar cycle, maximum amplitude, rising time, descending time 
PACS: 96.50.Wx, 96.60.-j, 96.60.Q-, 96.60.qd 
doi: 10.1360/SSPMA2015-00097 
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